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ABSTRACT 

Recent hydrostatic X-ray studies of the hot interstellar medium (ISM) in early-type 
galaxies underestimate the gravitating mass as compared to stellar dynamics, imply- 
ing modest, but significant deviations from exact hydrostatic equilibrium. We present 
a method for combining X-ray measurements and stellar dynamical constraints in the 
context of Bayesian statistics that allows the radial distribution of the implied non- 
thermal pressure or bulk motions in the hot ISM to be constrained. We demonstrate 
the accuracy of the method with hydrodynamical simulations tailored to produce a 
realistic galaxy model. Applying the method to the nearby elliptical galaxy NGC 4649, 
we find a significant but subdominant nonthermal pressure fraction (0.27 ± 0.06) in 
the central (< 5 kpc) part of the galaxy, similar to the level of deviations from hydro- 
static equilibrium expected in galaxy clusters. Plausible sources of systematic error, 
if important, may reduce this fraction. Our results imply ~360km s" 1 random turbu- 
lence or a magnetic field B = (39 ± 6)(n e /0.1 cm _3 )°- 59±0 09 /iG, whereas gas rotation 
alone is unlikely to explain the detailed nonthermal profile. Future observations with 
Astro-H will allow turbulence or gas rotation at this level to be detected. 



Key words: galaxies: elliptical and lenticular, cD — galaxies: 
netic fields — galaxies: kinematics and dynamics — turbulence- 



ISM — galaxies: mag- 
- X-rays: galaxies 



1 INTRODUCTION 

The distribution of mass in early-type galaxies is a funda- 
mental yardstick for understanding their formation and evo- 
lution. Our current ACDM cosmological paradigm predicts 
ubiquitous, massive dark matter halos, the distribution of 
matter within which correlates with the mass of the sys- 
tem (e.g. Navarro et al. 1997; Bullock et al. 2001; Maccio 
et al. 2008) . Within the optical radius of the galaxy, the rel- 
ative contributions to the gravitational potential of the dark 
and luminous matter may explain the tilt of the fundamen- 
tal plane (e.g. Cappellari et al. 2006; Bolton et al. 2007; 
Humphrey & Buote 2010), and provide insights into the 
poorly understood processes by which they interact gravita- 
tionally, in particular the interplay between adiabatic con- 
traction and dynamical friction (e.g. Blumenthal et al. 1986; 
Gnedin et al. 2004; El-Zant et al. 2004; Humphrey et al. 
2006; Dutton et al. 2007; Gnedin et al. 2007; Abadi et al. 
2010; Napolitano et al. 2010; Maccio et al. 2012). The stellar 
mass-to-light (M/L) ratio can provide important constraints 



on the history of star formation and the shape, and univer- 
sality, of the stellar initial mass function (IMF) (e.g. Ger- 
hard et al. 2001; Gnedin et al. 2007; Humphrey et al. 2009; 
Treu et al. 2010; Cappellari et al. 2012). Furthermore, at the 
smallest scales, the tight relations between the central super- 
massive black hole (SMBH) mass and the global properties 
of the host galaxy imply strong evolutionary symbiosis (e.g. 
Kormendy & Richstone 1995; Gebhardt et al. 2000; Merritt 
& Ferrarese 2001; Giiltekin et al. 2009; Ferrarese & Ford 
2005, for a review). 

A number of techniques are available to study the mass 
distributions in early- type galaxies, each having its own dis- 
tinct advantages and disadvantages. For example, gravita- 
tional lensing studies directly probe the projected mass dis- 
tribution, but alone they do not generally produce detailed 
mass profiles for individual giant elliptical galaxies, although 
they can provide important constraints on the profile aver- 
aged over multiple systems (e.g. Gavazzi et al. 2007), or 
when combined with stellar dynamics measurements (e.g. 



© 0000 RAS 



2 P. J. Humphrey et al. 



Koopmans et al. 2009; Barnabe et al. 2011). By combining 
information from multiple dynamical tracers, sophisticated 
(axisymmetric) orbit-based (stellar) dynamical models are 
now able to probe from within the sphere of influence of 
the SMBH out to tens of kpc (e.g. Romanowsky et al. 2003; 
Gebhardt & Thomas 2009), but models self-consistently in- 
corporating dark matter halos can be computationally costly 
to evaluate (e.g. Shen & Gebhardt 2010), and care must be 
taken to minimize systematic uncertainties associated with 
their implementation (e.g. Valluri et al. 2004; Shen & Geb- 
hardt 2010; Das et al. 2011; Long & Mao 2012) and the (gen- 
erally) unknown inclination of the system (Gavazzi 2005; 
Thomas et al. 2007). 

Hydrostatic X-ray techniques, on the other hand, are 
appealing in part due to their computational simplicity, 
given the isotropy of the gas pressure tensor and the fact that 
the spherical approximation (which is generally assumed, al- 
though not exclusively; e.g. Buote & Canizares 1994, 1996, 
1998; Buote et al. 2002; Statler & McNamara 2002) typically 
introduces only a small bias (Buote & Humphrey 2012c, 
and references therein), particularly if the spherically aver- 
aged mass profile is close to a singular isothermal sphere 
(Buote & Humphrey 2012b; Churazov et al. 2008), which is 
generally true in early-type galaxies (e.g. Koopmans et al. 
2009; Humphrey & Buote 2010; Churazov et al. 2010a). 
These methods are especially powerful since current instru- 
mentation allows measurements from close to the sphere of 
influence of the most massive nearby black holes out to 
^ R-2500 for some nearby galaxies (Humphrey et al. 2008, 
2009, 2011, 2012b; Wong et al. 2011; Buote & Humphrey 
2012a, for a review), and the same approach can be applied 
self-consistently from —Milky Way-mass galaxies to massive 
clusters. 

The accuracy of hydrostatic methods is, naturally, con- 
tingent on the extent to which equilibrium gas motions 
reflect the thermal gas pressure rather than turbulent or 
streaming motions. To ensure systems are not grossly far 
from equilibrium and do not exhibit dynamically dominant, 
supersonic gas motions, excluding objects with highly asym- 
metric X-ray isophotes is an important first step (Buote & 
Humphrey 2012a). Nevertheless, subsonic gas motions (ei- 
ther bulk or turbulent), or nonthermal support from mag- 
netic fields or cosmic rays could still be important (e.g. Chu- 
razov et al. 2008). Cosmic ray pressure is likely to be most 
important in the vicinity of cavities inflated by an AGN jet 
(e.g. Mathews & Brighenti 2008). Magnetic fields weaker 
than a few /xG in the centre of the galaxy should not be 
dynamically important (Buote & Humphrey 2012a), but 
there is considerable unccrtaintiy on the actual magnetic 
field strength, with current observational constraints from 
galaxies with embedded radio jets ranging from — 1/iG to a 
few tens of /iC (Vallee 2011). 

Interesting constraints on gas dynamics are few (Buote 
& Humphrey 2012a). The upper limits on X-ray spectral 
line broadening (~400-500km s" 1 , from XMM RGS obser- 
vations: Sanders et al. 2011), and velocity gradients (—1000- 
2000km s _1 from the spatial variation of the X-ray line cen- 
troids in clusters: Ota et al. 2007; Sato et al. 2008; Sugawara 
et al. 2009; Tamura et al. 2011) exceed the typical circular 
velocity for a giant elliptical galaxy. Indirect gas dynamical 
constraints in the very centres of early-type galaxies can be 
obtained from resonance scattering (or its absence) in strong 



emission lines (Churazov et al. 2010b, for a review). Cur- 
rent measurements suggest microturbulent velocities rang- 
ing from — 300-700km s _1 and 100-500km s _1 respectively 
in two giant elliptical galaxies (albeit sensitive to the fitted 
plasma code; de Plaa et al. 2012), to < 100km s _1 in the cen- 
tre of the elliptical galaxy NGC 4636 (Xu et al. 2002; Werner 
et al. 2009, who also found evidence of little turbulence in 
two other galaxies). The interpretation of these results, how- 
ever, may be complicated by the scale and anisotropy of the 
turbulence (Zhuravleva et al. 2011). Dynamically important 
gas rotation, such as expected in subsonically inflowing gas 
due to angular momentum conservation, should result in a 
flattening of the X-ray isophotes parallel to the rotation axis. 
This allowed Brighenti et al. (2009, hereafter B09) to infer 
— 200-300km s _1 gas motions in the central —kpc of the 
elliptical galaxy NGC 4649, and led them to suggest that 
inflowing, rotating gas may be common in the very central 
parts of nearby elliptical galaxies. 

A direct assessment of the accuracy of hydrostatic 
methods in an elliptical galaxy can be made by compar- 
ing the mass inferred from the X-rays to that obtained by 
independent techniques. Orbit-based stellar dynamical mod- 
els are the alternative method of choice for nearby systems, 
as they provide high-quality mass profiles that overlap well 
in radial scale with X-ray work. To date such comparisons 
have been made reliably in only a handful of cases. Geb- 
hardt & Thomas (2009) found inconsistency between their 
axisymmetric dynamical models for M87 and the gravita- 
tional potential inferred from X-rays by Churazov et al. 

(2008) . Das et al. (2010) and Murphy et al. (2011) quan- 
tified this discrepancy further, finding that the X-rays un- 
derestimate the optically inferred mass by —50% within the 
central —4 kpc, while the data generally agreed at larger 
scales. M 87 exhibits, however, substantial disturbances to 
its X-ray morphology (e.g. Bohringer et al. 1994; Forman 
et al. 2007; Million et al. 2010), which may introduce sys- 
tematic uncertainties into a single-phase spherical hydro- 
static mass analysis. Shen & Gebhardt (2010) found that 
the X-ray inferred mass that was measured by Humphrey 
et al. (2008) for the elliptical galaxy NGC 4649 was -40% 
lower than their dynamical models implied inside —8 kpc, 
despite its very smooth, relaxed X-ray isophotes. A similar 
result was obtained by Das et al. (2010), who also compared 
their X-ray determined mass profile for NGC4472 to (un- 
published) orbit-based models, finding it to be —30% lower 
inside —10 kpc but —20% higher further out. This dichotomy 
may, in part, reflect X-ray disturbances in the system at the 
largest scales (Irwin & Sarazin 1996; Biller et al. 2004). Rusli 
et al. (2011) measured the mass distribution in the central 
regions of the lenticular galaxy NGC 1332 with axisymmetric 
orbit-based models (excluding a dark matter halo), and find- 
ing the X-ray inferred profile obtained by Humphrey et al. 

(2009) to be -50% lower. 

A number of other studies are less definitive due to limi- 
tations in the optical or X-ray analysis. For example, Norris 
et al. (2012) derived the mass of NGC 3923 from axisym- 
metric orbit-based models, including a dark matter halo, 
finding an offset from the X-ray mass profile of Fukazawa 
et al. (2006). Although corrections were applied to account 
for unresolved sources when computing the temperature, 
Fukazawa et al. did not account for this contamination when 
deriving the density profile. An independent Rosat study by 
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Buote & Canizares (1998), considering non-spherical mass 
models for this galaxy, was actually in better agreement with 
the dynamical mass. Several authors have compared X-ray 
inferred masses with simpler dynamical models that may 
be less reliable (e.g. Mathews & Brighenti 2003b; Ciotti & 
Pellegrini 2004; Churazov et al. 2008, 2010a; Romanowsky 
et al. 2009; Johnson et al. 2009; Das et al. 2010), or have 
used X-ray mass profiles that are likely unreliable due to 
large morphological disturbances (e.g. NGC4636: Johnson 
et al. 2009; NGC5846: Churazov et al. 2010a; Das et al. 
2010), unresolved point source contamination (Pellegrini & 
Ciotti 2006), or systematic uncertainties associated with the 
model parameterizations (Romanowsky et al. 2009). 

Based on these studies, a consistent picture is begin- 
ning to emerge of X-ray hydrostatic methods underestimat- 
ing the mass inferred from stellar dynamical techniques by 
on average ~30-50% in the central part of the galaxy. Un- 
derstanding this discrepancy is essential not only in order to 
correct hydrostatic mass for this effect, but also to gain in- 
sight into the behaviour and fate of the gas at the centre of a 
galaxy-scale cooling flow (Mathews & Brighenti 2003a). Ob- 
taining detailed radial profiles of the measured nonthermal 
pressure or bulk gas motions is a critical next step. While, in 
essence, this can be achieved by differencing the mass pro- 
files inferred by both techniques, most recent studies have 
only compared them approximately (e.g. Shen & Gebhardt 
2010; Das et al. 2010, 2011; Rush et al. 2011; Murphy et al. 
2011; Norris et al. 2012) 1 

Given that the level of discrepancy is < 50%, it is cru- 
cial to control tightly the systematic uncertainties in both 
the X-ray and stellar dynamical measurements. From an X- 
ray perspective, there are a number of ways to infer the 
gravitating mass from X-ray data, as reviewed in Buote & 
Humphrey (2012a). For a reliable measurement, however, 
any fitted models must have sufficient flexibility to cap- 
ture the full range of physical gas temperature and density 
profiles. The traditional "smoothed inversion" technique, as 
widely employed in studying galaxy clusters, involves fit- 
ting parameterized models directly to the density and tem- 
perature distributions but, due to their diverse shapes in 
elliptical galaxies, there are in general no well-defined "uni- 
versal" parameterizations that can be employed, introducing 
potentially large systematic uncertainties in this regime (e.g. 
Humphrey et al. 2009; Buote & Humphrey 2012a). Das et al. 
(2010) proposed a minimally parametric "smoothed inver- 
sion" technique that does not a priori restrict the form of 
the density and temperature profiles, but it does employ a 
certain level of smoothing. We discuss the performance of 
this method in the Appendix. In practice, we have found 
an entropy-based "forward fitting" approach to be an at- 
tractive means for computing the mass profile reliably, as 
it rigorously enforces both a physical mass distribution and 
Schwarzschild's stability criterion for the ISM (Humphrey 

1 We note that if one compares the difference between two radii 
in the potential (rather than the enclosed mass) that is obtained 
with each method, as advocated by Churazov et al. (2008), the 
average nonthermal pressure in that range can also be derived. 
The interpretation is complicated, however, by the lack of statis- 
tical independence between the potential "data-points" generated 
in this way. Nevertheless, Churazov et al. used this approach to 
infer the global average level of support in two objects. 
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Figure 1. Comparison of the cicular velocity profile of NGC 4649 
derived from the stellar dynamical analysis of Shen & Gebhardt 
(2010, grey region), and our purely hydrostatic analysis of the 
hot ISM (blue shaded region). We also show (red) the best-fitting 
model and 1-cr confidence region obtained by Das et al. (2010) 
from their hydrostatic X-ray analysis. The modest discrepancies 
between the X-ray measurements likely arise from different mod- 
elling assumptions in the Das et al. analysis (see the Appendix), 
while the higher optical measurement can be interpreted as aris- 
ing from modest nonthermal pressure. 

et al. 2008, 2009, 2011). The drawback of this technique 
for computing the nonthermal pressure profile is that it as- 
sumes the mass inferred by hydrostatic techniques can be 
well-parameterized by the same model that can fit the true 
mass distribution, which need not be true in the case of a 
nonthermal pressure gradient. 

In this paper, we develop a simple extension of the 
entropy-based, forward fitting technique that incorporates 
non-hydrostatic effects. This method allows the radial dis- 
tribution of the nonthermal pressure or gas motions to be 
constrained directly from a joint, Bayesian analysis of the 
X-ray and stellar dynamics data. We illustrate the perfor- 
mance of the method both with simulated galaxies, and with 
a real system (NGC 4649), and interpret the results in terms 
of physical models for the gas flow, turbulence and possible 
magnetic fields. In what follows, all error bars correspond to 
1-cr uncertainty. 



2 NGC 4649 

To test our method, and provide important constraints on a 
well-studied system, we focused our analysis on NGC 4649. 
This is a nearby (15.7 Mpc) 2 , giant elliptical galaxy with a 
luminous X-ray halo and a remarkably round, relaxed X-ray 
morphology, at least within ~20 kpc (Trinchieri et al. 1997; 
Randall et al. 2004; Humphrey et al. 2008). Low significance 
surface brightness asymmetries were reported in a shallow 

2 We adopt a distance of 15.7 Mpc for consistency with the dy- 
namical models of Shen & Gebhardt (2010) 
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Chandra observation by Randall et al. (2004) , which they, in 
part, speculated may arise due to convective flows. This hy- 
pothesis is inconsistent, however, with the monotonically ris- 
ing entropy profile (implying global stability against convec- 
tion) . In any case, the features were not confirmed in deeper 
Chandra observations (Humphrey et al. 2008). Shurkin et al. 
(2008) reported the detection of small cavities associated 
with the weak radio jet but these were also not confirmed 
in our deeper data. Instead the gas was found to be very 
symmetrical and relaxed-looking (Humphrey et al. 2008). 
The X-ray isophotes do exhibit central flattening, however, 
which suggests that gas rotation may be important in the 
central ~1 kpc (B09). 

Given its proximity and high surface brightness, we were 
able to measure the gas properties from close to the sphere 
of influence of the central black hole (~1") out to ~25 kpc, 
allowing, for the first time, the mass of any SMBH to be in- 
ferred directly from a hydrostatic X-ray model (Humphrey 
et al. 2008). State of the art, orbit-based, axisymmetric stel- 
lar dynamical models, which spanned a similar radial range, 
were constructed by Shen & Gebhardt (2010), assuming an 
edge-on geometry (but verifying that the results were not 
very sensitive to the adopted inclination). Although the in- 
ferred black hole mass was in statistical agreement with the 
X-ray measurement, they found a modest discrepancy with 
the hydrostatic mass profile within the central ~8 kpc. A 
similar discrepancy was noted by Das et al. (2010) between 
the dynamical mass of Shen & Gebhardt and their own X- 
ray mass model. 

In Fig 1 we show the circular velocity (v c = \J GM/r, 
where G is the Universal gravitational constant and M is the 
total mass enclosed within radius r) profile for NGC 4649, 
based on the X-ray data (Humphrey et al. 2008) and the 
dynamical analysis of Shen & Gebhardt (2010) 3 , illustrat- 
ing the offset between the inferred mass distributions. We 
also overlay the v c profile of Das et al. (2010), which overall 
agrees well with our measurement over the range ~3-9 kpc, 
but is marginally higher outside this range. We ascribe these 
differences to particular choices made in their analysis (see 
the Appendix for a more detailed discussion). The generally 
good agreement between the two X-ray measurements indi- 
cates that the level of discrepancy with the optical results is 
not sensitive to the details of the X-ray modelling although, 
as discussed in the Appendix, the error bars may be. 

In the present study, we interpret the offset between the 
X-ray and dynamical measurements as wholly arising due 
to deviations from hydrostatic equilibrium. Naturally, this 
is contingent upon the optical result being an unbiased mea- 
surement of the true gravitational potential. For the state of 
the art models used by Shen & Gebhardt (2010), which have 
previously shown to perform well in recovering the mass of 
simulated galaxies, we do not expect there to be significant 
inherent systematic uncertainties, other than a dependence 
on the viewing angle (Thomas et al. 2007). We therefore 
assume that the gravitational potential inferred by Shen & 
Gebhardt (2010) is accurate. We note that, recently, Das 
et al. (2011) constructed "made to measure" particle-based 



3 This Bayesian realization of the v c profile was derived by lin- 
early interpolating over the x 2 -topology obtained by Shen & Geb- 
hardt (2010) from their stellar dynamical analysis. 



models for the system (using the same kinematics data in 
the central region), finding larger statistical errors than Shen 
& Gebhardt, and thus formal consistency with the X-ray 
profile. However, investigating the origin of such systematic 
differences is beyond the scope of the present paper. 

In order to constrain the nonthermal pressure profile, we 
carried out a joint analysis of the optical and X-ray data, 
as described in § 3.2. For the X-ray data, we used the pub- 
lished gas density and temperature profiles from Humphrey 
et al. (2008). For the optical constraints, we used the poste- 
rior probability distribution obtained by Shen & Gebhardt 
(2010) and linearly interpolated over their grid of x 2 values 
(and assumed flat priors for each parameter). 



3 MEASURING NONTHERMAL PRESSURE 
3.1 Preliminaries 

If the ISM is in an equilibrium state and the magnetic fields 
are tangled, the total gravitating mass, M, enclosed within 
a surface S is given by (Fang et al. 2009): 



M = 



1 

4^G 



Pg 



(v.V)v 



.dS 



(1) 



where G is the universal gravitational constant, p g is the 
gas density, P is the total pressure (including thermal gas 
pressure, and pressure from tangled magnetic fields, cosmic 
rays and turbulence), v is the (time averaged) gas velocity 
field and dS a surface element. Assuming S is a spherical 
surface of radius r, this becomes: 



Vc_ 

r 



1 

47T 



J_dP 

Pg dr 



dQ 



4t£ 



f. (v.V)v dQ, 



(2) 



where v c is the circular velocity, Q is the solid angle, and f 
is the unit vector in the radial direction. For convenience, 
we define < p g > such that 



d<P> 



dr 



i_dp 

Pg dr 



dQ 



(3) 



< Pg > » «• j 4n 

where < P > is the spherically averaged pressure. We also 
define an "effective" gas velocity, v e f / , such that 

2 

f. (v.V)vdft (4) 



v: 



ff 



1 

47T 



477 



To make further progress, we assume that the density 
and pressure are approximately stratified on concentric 
spheres, which is generally a good approximation (Buote & 
Humphrey 2012b, c, and references therein). In that case, we 
can drop the <> averaging notation in Eqn 3. If the only dy- 
namically important gas motions are rotation about a com- 
mon axis (i.e. v = v^e^ where e^ is the unit vector in the 
4> direction in spherical coordinates), v\jj = j^v^dQ/Aiv; 
hence v e f / would be a spherically averaged gas rotation ve- 
locity. 

We define P = P a /(1 — fnth), where P g is the ther- 
mal gas pressure (P g — p g kT '/ '(jimj); k is Boltzmann's 
constant, T is the gas temperature, /i ~ 0.62 is the mean 
molecular weight, and tuh is the weight of hydrogren), and 
fnth is the nonthermal pressure fraction. We also define 
S — p g 2 ^ 3 kT/(/j,mH), which is proportional to the tradi- 
tional definition of the entropy proxy, K — n^ 2 ^ z kT, where 
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n e is the electron number density. Combining Eqns 2-4 and 
rearranging, we obtain 



S 



d_ 

dr 



1 — frith 



= -r [v 2 c 



2 

V *ff 



Writing vt = v1. g + v'f 



c>ng , where v c , g 



(5) 

is the circular veloc- 



ity due to the gas mass, and v c ,ng is due to the non-gas 
mass, differentiating and folding in the gas mass continuity 
equation, 



J- (rvl g ) = AtvGt 2 
we obtain: 



(6) 



d_ 

dr 



s 



d_ 

dr 



1 - fn 



dr 



-4nGr 



V eff] 



V 5 



(7) 



3.2 Method 



Adopting parameterized models for S, v Ct , lg , f n th and v e ff, 
Eqn 7 can be integrated numerically, to solve for P using a 
fourth-order Runge-Kutta method. We start the integration 
at some arbitrary, small radius ro. Two boundary condi- 
tions are needed; we set the thermal gas pressure at ro to 
a value to be determined by our fit, and assume that the 
gas mass is zero within ro. Since the model profiles of P g 
and S uniquely determine the gas density and temperature, 
the observed temperature and density profiles can be fitted 
by adjusting the various model parameters. As the temper- 
ature and density are generally inferred by spherical depro- 
jection of spectra accumulated in a series of concentric bins, 
to ensure a reliable comparison we integrated appropriately 
weighted functions of pressure and temperature to compute 
the predicted contribution to the emission measure (and 
hence mean density) from gas in the corresponding shell, 
as well as an emission-weighted temperature. We compared 
these averaged quantities to our measured temperature and 
density data points, as discussed in Humphrey et al. (2008) 
and Humphrey et al. (2009). See Appendix B of (Gastaldello 
et al. 2007b) for a detailed discussion of incorporating the 
plasma emissivity into the gas modelling. 

To parameterize v Ct , lg , we model the gravitating mass 
as arising from a central black hole, a stellar mass com- 
ponent (assuming mass follows light), plus a dark matter 
halo, which can either be a Navarro et al. (1997, hereafter 
NFW) profile, or a "cored logarithmic" potential (Binney 
& Tremaine 2008). The black hole mass, stellar M/L ratio, 
and normalization and characteristic scale of the dark mat- 
ter component are adjustable fit parameters. To parameter- 
ize S, we adopt a multiply broken powerlaw, plus a constant 
term. This model has sufficient flexiblity to capture the over- 
all shapes of the entropy profiles in galaxy groups and clus- 
ters (e.g. Gastaldello et al. 2007a; Sun et al. 2009; Cavagnolo 
et al. 2009; Humphrey et al. 2012a). Since a priori there is 
no obvious parameterized form for f nt h or v e ff, we parame- 
terize them with a cubic spline function. While the number 
of spline knots is arbitrary, we found the model to be gener- 
ally well-behaved if there are roughly half as many knots as 
there are density (or temperature) bins, and they are placed 




R (kpc) 

Figure 2. Recovered v e ee for the hydrostatic simulation. The 
data-points correspond to the spline knots fitted to the simulated 
data, and the grey region indicates the 1-cr confidence region of 
the model interpolated between them. The magenta line is the 
circular velocity (v c ; set to match the X-ray mass inferred by 
Humphrey et al. 2008), and the arrows indicate the approximate 
centre of each annulus used in the spectral analysis. The true 
v e ff distribution (v e ff= 0) is consistent with the best fitting 
model, although there is a slight (~30km s _1 ) bias, which is not 
dynamically important. 



approximately logarithmically in radius. Since f n th only ap- 
pears inside the derivative in Eqn 7, a boundary condition 
on it is required. We therefore set f n th= at large radii. 

In the hydrostatic approximation, / = and v e ff = 
and so the two unknowns (S and v c ) can be constrained 
completely by the two observables (p g and kT). By explor- 
ing parameter space with a Bayesian code (specifically, we 
use version 2.7 of the MultiNest code 4 (Feroz et al. 2009)), 
we can also fold in stellar dynamical constraints on v c by us- 
ing the posterior probability distribution from the dynamical 
analysis as priors on some of the parameters. In this case, 
we can relax the hydrostatic approximation and addition- 
ally measure the profile of either f nt h or v e ff. These two 
quantities are highly degenerate and so, without additional 
constraints, they cannot be disentangled uniquely. We there- 
fore adopt the pragmatic approach of fitting each profile in 
turn, while setting the other term to zero. 



3.3 Tests 

To test the ability of our method to recover the nonthermal 
pressure profile in early- type galaxies, we used it to investi- 
gate two simulated datasets, each having different amounts 
of nonthermal support. Since our method involves directly 
fitting the X-ray data, and since we assume that the stel- 
lar dynamical constraints provide an unbiased estimate of 
the true gravitational potential, we considered the limit of 
perfect optical constraints, but a typical X-ray exposure 
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Figure 3. Left: Marginalized v e ff from fitting the B09 simulation. The data-points indicate the spline knots fitted to the simulation, 
and the grey region is the 1-cr confidence region of the model interpolated between them. The blue line indicate the true profile of v e ff. 
The magenta line is the circular velocity (matching the mass inferred by Humphrey et al. 2008), and the arrows indicate the approximate 
centre of each annulus used in the spectral analysis. Overall, we find good agreement between the true distribution and the model. Right: 
The same nonthermal support profile parameterized as a nonthermal pressure fraction (f n th \ g rev shaded region and data-points), shown 
with the true profile (blue) and (arrows) the approximate centres of each radial bin. The overall fit is also fairly good . We note that 
the modest discrepancy between the true and inferred f e // profiles around ~0. 1-0.2 kpc occurs between two spline knots and is mostly 
confined to the central radial bin, so we caution against over-interpretation. 



(100 ks). The first of these datasets was derived assum- 
ing spherical symmetry and perfect hydrostatic equilibrium 
(i.e. v e ff=f nt h= 0; hereafter the "hydrostatic simulation"). 
The second simulated dataset was based on the axisymmet- 
ric, hydrodynamical simulation performed by B09, in which 
there is turbulence and bulk gas motion. In both cases, the 
simulations were tailored to match approximately the prop- 
erties of NGC 4649, as observed with Chandra. 

For each simulated dataset, we started with the tem- 
perature, density and abundance of the g function 
of position in an axisymmetric (R,z) grid, and created an 
artificial Chandra events file using a Monte Carlo method 
to generate photons. We assumed the source is at 15.7 Mpc, 
and the observation was 100 ks in duration. At each grid po- 
sition, a random number of photons were generated, assum- 
ing an APEC plasma with the appropriate temperature and 
gas abundance, and using representative Chandra responses 
for close to the aimpoint. Each photon was randomly as- 
signed a PHA (energy) bin and a spatial position within the 
bin, before being projected onto the sky. We processed the 
simulated events file similarly to the real Chandra data of 
NGC 4649 (Humphrey et al. 2008), which involves spherical 
deprojection with the pro jet model in Xspec vers 11, and fit- 
ting an APEC plasma model with variable abundances, and 
produces the temperature, density and abundance profiles. 
Where appropriate, we tied the abundance between adjacent 
annuli, to improve constraints. To simplify the analysis, we 
filtered out photons originating from outside the outermost 
shell used in our deprojection. 

We fitted the resulting temperature and density profiles 
simultaneously with our modified hydrostatic model, assum- 
ing flat priors on all fit parameters. In these tests, we aimed 
to determine how accurately the nonthermal pressure profile 
can be recovered in the limit of perfect dynamical data, so we 



fixed the mass model parameters to match the true distribu- 
tion used in the simulations. For the hydrostatic simulation, 
we show the recovered distribution of v e fj in Fig 2, demon- 
strating that the true distribution (v e ff= 0) was reasonably 
well recovered. There was a slight (~ 30km s _1 ) bias in the 
best-fitting velocity, which may reflect slight uncertainties in 
the deprojection procedure, but such a small velocity would 
be dynamically unimportant and so has minimal effect on 
our conclusions. 

In Fig 3 (left panel) we show the fitted v e ff profile for 
the B09 simulation, along with its true distribution. Since 
the velocity field predicted by the simulation was known, we 
were able to compute the true v e f / profile by explicitly eval- 
uating Eqn 4 (for more details of how this was done, see Fang 
et al. 2009). The fitted model was, by definition, a smooth 
function that cannot exactly capture all of the features in 
the true profile, but the overall shape was well recovered, 
especially in the vicinity of each spline knot. The largest de- 
viations actually occurred within the central radial bin, at 
scales that are, effectively, unresolved. In the inner regions of 
the simulation (< 1 kpc) the X-ray isophotes were quite flat- 
tened and there was a strong azimuthal temperature gradi- 
ent (B09), so the spherical, single-phase approximation that 
is implicit in the deprojection method we adopted is a simpli- 
fication. Despite this, the impact of these effects was modest 
and the recovered w e // profile was sufficiently accurate for 
our purposes. 

In Fig 3 (right panel), we show the fitted fnth profile 
for the B09 simulation, and its true distribution. The non- 
thermal support in the simulation actually arises (mostly) 
from bulk gas motions rather than an additional pressure 
component. Still, f nt h and v e // are highly degenerate, so we 
can still parameterize the nonthermal support as an "effec- 
tive" f„th- Writing out Eqn 5 for the two cases (« e //= 
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Figure 4. Profiles of kT (upper panel) and gas density (lower 
panel), along with the best-fitting model. Note the good agree- 
ment in all radial bins. 



and f n th= 0), and eliminating dP g /dr between them, we 
find that f„th and v e ff are related by: 



2 

V eff 



. 2 kT dlog(l -/„«,) 

JnthV c + 



(8) 

fimn d log r 

We integrated this equation numerically to obtain the "true" 
fnth , which we found to be in good agreement with the fitted 
profile (Fig 3). 



4 APPLICATION TO NGC 4649 
4.1 Results 

To study NGC 4649, we adopted the same parameterized 
gravitating mass model as Shen & Gebhardt (2010) (black 
hole, stellar component and a cored logarithmic dark matter 
mass model) , and folded in the results of the stellar dynami- 
cal analysis as priors on the black hole mass, stellar M/L ra- 
tio and parameters of the dark matter model. For the stellar 
component, Shen & Gebhardt deprojected the V-band stel- 
lar light distribution from Kormendy et al. (2009), assum- 
ing an edge-on, oblate spheroidal geometry and a constant 
axis ratio of 0.9. Unlike the axisymmetric assumption of the 
dynamical analysis, the X-ray technique assumes spherical 



geometry, and so it is important to average this profile spher- 
ically (Buote & Humphrey 2012b, c), which we did numeri- 
cally. 

We separately parameterized v e ff and f n th as a 7- knot 
spline (each knot being placed approximately logarithmi- 
cally in radius), and assuming flat priors on the model val- 
ues at each knot point, within reasonable limits. The best- 
fitting model actually captured the shapes of the density 
and temperature profiles acceptably, although the f„th pa- 
rameterization is preferred (x 2 =18.7/13 and 22.1/13 for the 
fnth and v e ff paramctcrizations, respectively; see Fig 4). 
The resulting v e ff and f n th profiles are shown in Fig 5. 

The shape of the v e ff profile is relatively flat, and 
could be consistent with a constant gas rotation velocity 
( v eff =232 ± 37km s _1 ), or with a mildly falling profile that 
becomes consistent with zero outside ~5 kpc. Expressed in 
terms of fnth, the profile is also consistent both with a con- 
stant nonthermal pressure fraction (f n th= 0.27 ± 0.06) or 
with a gradually declining distribution that becomes con- 
sistent with zero outside ~2 kpc. Still, within ~2 kpc, ap- 
proximately ~30% of the support must come from bulk gas 
motions or nonthermal pressure. 

4.2 Systematic error budget 

All methods for inferring the mass distributions in early- 
type galaxies entail arbitrary analysis choices, some of which 
may quantitatively affect the measured mass. In this section, 
we briefly explore how sensitive ourn conclusions are to the 
choices in the X-ray modelling. Further systematic uncer- 
tainties may be associated with choices in the stellar dynam- 
ical analysis (for example, the inclination angle assumed), 
but exploring them is beyond the scope of the present pa- 
per. The systematic error budget is summarized in Table 1, 
and we discuss below how each test was performed. For a 
more detailed discussion of the various systematic error as- 
sessments considered here, see Humphrey et al. (2012a). 

We first examined the sensitivity of our spectral- 
fitting results to the treatment of the Chandra background 
("ABackground") by using the standard "background tem- 
plate" spectra, suitably renormalized to match the data 
at > 10 keV, instead of the more robust, modelled back- 
ground adopted by default (Humphrey et al. 2008). We 
experimented with replacing the APEC thermal plasma 
model used in spectral fitting with the MEKAL model 
("ASpectral") and varying the distance by ^33% (corre- 
sponding to the statistical error on the distance to NGC 4649 
given by Tonry et al. 2001, "ADistance"). To investigate how 
changing the parameterized model used to fit the v e ff or 
fnth profiles can affect the results, we experimented with re- 
ducing the number of spline knots to 3 (at 0.5 kpc, 0.16 kpc 
and 20 kpc), and also with using a constant v e ff or f n th 
("AKnots"). These simpler parameterizations did not lead 
to significantly poorer fits. Next we examined how the choice 
of priors might be influencing our results by replacing each 
of the flat priors we adopted in our analysis by priors which 
were flat in logarithmic space ("APriors"). We also con- 
sidered changing the entropy parameterization by adding 
an additional break at large radii ("AEntropy"), experi- 
mented with folding in the full covariance between adjacent 
temperature and density data-points ("ACovariance") and 
turning off the emissivity computation when evaluating the 
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Figure 5. The same as Fig 3 but for the real NGC4649 data. In this case the "true" profile is not shown as it is a priori unknown, and 
v c is derived from the joint fit to the stellar dynamical and X-ray data. 




R (kpc) n e (cm 3 ) 

Figure 6. Left panel: The implied magnetic field profile in NGC4649, assuming the nonthermal pressure comes from the B-field in 
equipartition with ~230 km s" 1 random turbulence. The data-points are the fitted spline points, while the grey region is the 1-a 
confidence range, interpolating between them. The arrows indicate the centres of each temperature and density data-bin. Right panel: 
The corresponding B-field versus electron density. The solid red line is the best- fitting powerlaw relation (see text). The dashed blue line 
is the prediction for the model of Mathews & Brighenti (1997, their Fig 2). The model slope agrees well with the observed data. The 
quoted uncertainties do not include systematic errors, including those from the stellar dynamical analysis. In particular, we note that, if 
the dynamical analysis of Das et al. (2011) is taken at face value, the implied level of nonthermal support may be lower than inferred 
here (see text). 



weighted temperature and density models ( "A Weighting" ; 
see Humphrey et al. 2012a for more details). We list the 
largest changes and statistical errors associated with each 
test in Table 1. In almost all cases, the change in the 
marginalized value was not much larger than the statistical 
error; in a few cases (especially at ~20 kpc) the statistical 
error-bars were enlarged by the test. Ultimately, none of the 
statistical uncertainties were large enough to lead to quali- 
tatively different conclusions. Nonthermal support, approx- 
imately at the ~25% level, is needed to reconcile the X-ray 
data with the published dynamical mass of NGC4649. 



5 DISCUSSION 

5.1 A non- hydrostatic X-ray model 

We have demonstrated that a simple modification to 
the entropy-based hydrostatic, forward fitting technique 
(Humphrey et al. 2008, 2009) allows nonthermal pressure 
profiles to be inferred from X-ray data, provided additional, 
unbiased, constraints on the mass profile are available (e.g. 
from stellar dynamical mass measurements). At scales re- 
solved by the density and temperature profiles, the new 
technique was able to recover reasonably well the radial dis- 
tributions of f n th or v e ff in a simulated galaxy, which was 
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Test 


v e ff(0.5 kpc) 
(km s~ 1 ) 


v ef f(1.6 kpc) 
(km s -1 ) 


v eff (20 kpc) 
(km s^ 1 ) 


fnth (0.5 kpc) 


/ nt h(1.6 kpc) 


/„ th (20 kpc) 


Marginalized 
Best-fit 


191 ± 57.3 
(230) 


(259) 


125 +77.8 

(45.9) 


n OO+0.06 
U-28_ 08 

(0.26) 


0.31 ± 0.06 
(0.31) 


00+0.15 
u - uu -o.ooo 
(0.00) 



ABackground 

ASpectral 

ADistancc 

AKnots 

AFit priors 

AEntropy 

ACovariancc 

AWeighting 



-34.6 



+25.8 (±63) 

/+90.6 
1,-24.9 
+20.2 /+53.81 
-67.9 1,-107 ) 

+58.2 (till 

+33.3 



2Z5 (i 93 - 1 ) 
+1.69 (±60.6) 

9 1+99. 9\ 
5.1 V-59.3J 

9.86 (±49.4) 



-11.4 (±55.3) 

+18.6 (±<°;|) 

-97.9 (±54.7) 
+ 11.7 /+41 
-20.2 V — 53.E 
+ 126 /+47.E 
-12.6 1,-112 
58.2 



+29 
-35 



-10.8 
±19.2 
+47.8 




0.025 (±0.06) 

0.05 



-26.4/ 
(+77.l\ 

-54.51 
f+36.8\ 

-53.9 



-67 (tH%) 
-63.7 (HZ ) 



+0.03 
+0.03 
+0.02 

+0.03 
-0.02 

+0.01 

+0.03 
-0.00 

-0.013 



■0.09 
-0.05 
■0.08 
-0.05 
0.09 
-0.08 
-0.11 
-0.04 

V -0.10 

/+0.09 
1-0.12 
/+0.06 
1,-0.07 



-0.077 (±0.06) 
+0.02 (±0.07) 
t°;°° (±0.07) 

-0.043 (±0.07) 

±8:SS ( ±tu °) 

+0.03 (±0.09) 



r0.05 



+0.02 (t° 



/+0.10 

0.07 
06 
07 



1-0. 18\ 
v-O.OoJ 
f +0.13\ 
-0.00 J 
/+0.29\ 
-0.10 J 



±0 

+0.003 

+0.10 

+0.28 (±0.07) 

+0.06 
+0.006 
+0.001 
±0 



1-0.06 ) 
/+0.16\ 
I— 0.01 J 
/+0.05N 

1-o.ooJ 

/+0.14\ 
l-O.OOi 



Table 1. Marginalized values and 1-cr confidence regions for v e ff and fnth measured at three representative radii. Since the best- fitting 
parameters need not be identical to the marginalized values, we also list the best-fitting values for each parameter (in parentheses). 
In addition to the statistical errors, we also show estimates of the error budget from possible sources of systematic uncertainty. We 
consider a range of different systematic effects, which are described in detail in § 4.2; specifically we evaluate the effect of treatment 
of the background (ABackground), spectral- fitting choices (ASpectral), distance uncertainties (ADistance), the number of spline knots 
used in the calculation (AKnots), priors on the fit parameters (AFit priors), the entropy parameterization (AEntropy), treatment of 
the covariance between data-points (ACovariance) and disabling the emissivity computation (AWeighting). We list the change in the 
marginalized value of each parameter for every test and, in parentheses, the statistical uncertainty on the parameter determined from 
the test. Note that the systematic error estimates should not in general be added in quadrature with the statistical error 



tailored to match real observations of the representative, 
bright system NGC 4649. At unresolved scales, the recovered 
distributions of v e ff and f n th were approximately correct, 
but were sensitive to the exact parameterizations used to fit 
them. 

The approach outlined here did not allow us to break 
the degeneracy between bulk gas motions (v e ff) and non- 
thermal pressure (fnth), and so we advocated determining 
each profile separately. In a real system, it is likely that 
some combination of both effects could contribute to non- 
thermal support, and so it is desirable to fold in additional 
constraints to try to disentangle them. This can be easily 
implemented in a Bayesian framework by modifying the pri- 
ors appropriately. In principle, resonance scattering mea- 
surements, or spectral line broadening constraints, could be 
employed in this way as restrictions on the turbulent pres- 
sure component of fnth- Similarly, magnetic field rotation 
measure constraints could be employed to restrict the mag- 
netic pressure contribution to f n th- As discussed in B09, 
the cllipticity profile of the gas actually provides impor- 
tant information on the gas motion, potentially allowing it 
to be used as an indirect constraint on v e ff- To do this 
fully self-consistently, as in B09, involves running a suite of 
hydrodynamical models, making it impractical to cover a 
suitably large region of parameter space. With significant 
simplifications, for example assuming that the only gas mo- 
tions are rotational, the problem becomes tractible. Unfortu- 
nately, to evaluate the models properly then involves relax- 
ing the spherical approximation, which substantially com- 
plicates the analysis beyond the scope of this paper. Non- 
spherical hydrostatic models including rotation have been 
constructed before (e.g. Buote k, Canizares 1996; Statler & 
McNamara 2002), but not for computing detailed tempera- 
ture and density distributions. 

As presented here, the model provides a means for mea- 
suring deviations from hydrostatic equilibrium in the gas. 



Conversely, if a physical model exists for v e ff or f n th, the 
same approach could be used to correct the X-ray mass mea- 
surement for deviations from hydrostatic equilibrium, mak- 
ing X-ray studies alone of nearby galaxies potentially suit- 
able for high-precision cosmology. At the moment, our un- 
derstanding of nonthermal support in the centres of galaxies 
is still in its infancy, but when the v e ff and f n th profiles of 
sufficiently large a sample of systems have been measured, 
this may become routine. 



5.2 Nonthermal support in NGC 4649 

Applying our model to NGC 4649, we have carried out a 
joint stellar dynamical and X-ray mass analysis, effectively 
reconciling the measurements made by Shen & Gebhardt 
(2010) and Humphrey et al. (2008) with a modest but sig- 
nificant (~25%) nonthermal pressure profile. We obtained 
these estimates by assuming that the ISM is maintained in 
static equilibrium, whereas in practice, we expect there to 
be a modest cooling flow (e.g. B09). The measured level 
of nonthermal support will then be an overestimate, albeit 
only very slightly if the cooling flow is very subsonic. By 
comparing the mass profiles obtained from stellar dynami- 
cal and X-ray analysis (Fig 1), Das et al. (2010) similarly 
concluded that the nonthermal pressure in NGC 4649 is no 
larger than ~32%, although they did not present a detailed 
fnth profile. This level of deviation from hydrostatic equilib- 
rium is fully comparable to that expected in galaxy clusters, 
which are routinely used for cosmology (e.g. Tsai et al. 1994; 
Buote & Tsai 1995; Evrard et al. 1996; Rasia et al. 2006; Na- 
gai et al. 2007; Piffaretti & Valdarnini 2008; Fang et al. 2009; 
Lau et al. 2011), and is consistent with (albeit slightly larger 
than inferred in) the conclusions of Churazov et al. (2008) 
for two more morphologically disturbed systems. Consider- 
ing the level of statistical and systematic errors in our pre- 
vious studies of early-type galaxies (e.g. Humphrey et al. 
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2006, 2008, 2009, 2011, 2012a,b), this level of nonthermal 
pressure, if ubiquitous, would not imply qualitatively differ- 
ent conclusions on the gravitating mass profiles. 

Since the level of nonthermal support is inferred from 
the subtle differences between the mass profiles derived with 
two techniques, control of the systematic errors is important 
for a robust measurement. Certain choices, such as the treat- 
ment of the X-ray background, the adopted fit priors, and 
the number of spline knots used, can have a non-negligible 
impact on the recovered profile, as is clear from § 4.2. Given 
present uncertainties (both statistical and systematic), the 
data currently do not allow us to determine whether the non- 
thermal support profiles (v e ff or fnth) are constant with ra- 
dius (as assumed by Churazov et al. 2008), or if they are 
gradually declining; they are clearly consistent with zero 
outside ~5 kpc. While we anticipate improved statistical X- 
ray constraints if we fold in the available archival XMM data 
and the significantly deeper Chandra observations (~200 ks, 
as compared to the 81 ks used in our analysis) that will soon 
be publicly available in the archive, the large statistical un- 
certainty on the current optically inferred mass at large radii 
(Fig 1) is a major limiting factor. 

At present, we have not folded in any systematic errors 
in the stellar dynamical analysis into our study. In particu- 
lar, uncertainties in the inclination of the system may be im- 
portant (e.g. Thomas et al. 2007). Shen & Gebhardt (2010) 
explored different inclination angles with a simplified model 
that included a fixed dark matter halo, finding that inclina- 
tion had little impact on their results. Given the restricted 
parameter space explored in this way, we did not attempt 
to incorporate these fits onto our work and marginalize over 
inclination. More intriguing are the recent particle based 
models of Das et al. (2011), who found much larger con- 
fidence regions for the enclosed mass profiles than Shen & 
Gebhardt when fitting kinematic data for the stars and plan- 
etary nebulae, implying a lower level of nonthermal support 
than measured here. Some of the discrepancies with the Shen 
& Gebhardt results may arise in part due to there being dif- 
ferent equilibrium configurations for the globular cluster and 
the planetary nebulae populations. However, this should not 
be important in the central part of the galaxy, where the in- 
ferred nonthermal pressure is strongest. This is supported by 
the good agreement between the stellar M/L ratios inferred 
by Shen & Gebhardt and Gebhardt et al. (2003), who used 
only stellar kinematics. While we expect the state of the art 
dynamical modelling of Shen & Gebhardt to be robust, it 
is unclear why the code used by Das et al. produced such 
a different result. Until this tension is resolved, some ques- 
tions persist over the exact level of nonthermal support in 
NGC 4649. As it stands, the true level of nonthermal support 
may be lower than we infer. Nevertheless, in the following 
sections, we consider the physical implications of nonther- 
mal support at the level inferred in our study. 

5.2.1 Gas rotation 

Given the smooth, relaxed isophotes of NGC 4649, we would 
not expect large-scale bulk gas motions other than, perhaps, 
gas rotation, which is expected in the centre of a cooling flow 
(e.g. Nulsen et al. 1984; Kley & Mathews 1995; Brighenti 
& Mathews 1996). B09 computed detailed hydrodynamical 
models for the gas flow in NGC 4649, which they tailored to 



match the gas density, temperature and ellipticity profiles. 
They concluded that rotation could be dynamically impor- 
tant within the central ~kpc. Comparison of their inferred 
v e ff profile (Fig 3, left panel) to the real data (Fig 5, left 
panel) strongly suggests that, while rotation could domi- 
nate the nonthermal support at these small scales, it cannot 
explain the observations at larger scales (~l-3 kpc). It is 
likely that B09 underestimated the level of gas rotation, as 
they adopted the gravitational potential inferred from the 
hydrostatic fit, rather than the deeper (true?) potential from 
dynamics, but the modest underestimate in v e ff that is im- 
plied in order to maintain the same level of isophotal flat- 
tening (< 20%) will still not reconcile the simulated rotation 
velocity with the observed profile at these scales. 



5.2.2 Random turbulence 

Alternatively, the nonthermal pressure profile {fnth) could 
arise due to turbulent motions within the gas. Assum- 
ing isotropic turbulence and eddy velocities drawn from 
a Maxwell-Boltzmann distribution, the turbulent pressure 
is given by P t ~ (l/3)p g Vt urb , where v tur b is the r.m.s. 
turbulent vel ocity. Adopting kT=0.8 keV, we find that 
vturb ^ 600^// n th/(l - f nth ) km s" 1 , or ~360±60 km s" 1 
for f n th~0.27 ± 0.06. Random turbulence of this magnitude 
may help to generate a dynamically important magnetic 
field that would also contribute nonthermal pressure (see 
§ 5.2.4). Whether plausible sources of random turbulence, 
such as mixing from stellar mass loss and Type la super- 
novae, galaxy merging and "sloshing", or AGN-driven dis- 
turbances in the ISM are sufficient to maintain turbulence 
at this level is unclear. 

Turbulent diffusion can provide a means for transport- 
ing angular momentum outwards, counteracting the ISM ro- 
tation that is expected to be induced by inflow, and thus 
making the X-ray image rounder (Brighenti & Mathews 
2000). Although sensitive to the assumed length-scale of 
the largest turbulent eddies, B09 used the isophotal flat- 
tening in the core of NGC 4649 to infer turbulent velocities 
~50 km s _1 , which are not dynamically important. 

It is worthwhile to compare the inferred level of turbu- 
lence with that implied in the centres of elliptical galaxies 
from other means. Sanders et al. (2011) placed an upper 
limit of 700km s _1 on line broadening in the XMM RGS 
spectrum of NGC 4649, which is too large to be interest- 
ing. While Werner et al. (2009) were not able to place use- 
ful constraints on resonance scattering in NGC 4649, they 
measured the effect in three other similar galaxies, includ- 
ing NGC 4636, for which they inferred v tur b < 100km s _1 . 
Conversely de Plaa et al. (2012) inferred significant turbu- 
lence (~140-700km s~ ) in two other systems, which may be 
consistent with our observation of NGC 4649. However, we 
note that, for one of these systems, NGC 5044, Sanders et al. 
(2011) concluded that the width of the spectral lines in RGS 
spectra was fully consistent with the source just being spa- 
tially extended, leaving little room for turbulent broadening. 
Whether turbulence as large as ~360km s _1 is realistic in 
the centre of a giant elliptical galaxy therefore remains un- 
clear at this point, but it does not appear ubiquitous. Still, 
this may be resolved by observations with Astro-H (§ 5.3). 
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5.2.3 Cosmic Ray Pressure 

Another possible cause of nonzero f nt h is cosmic ray pres- 
sure. Strong cosmic ray injection into the ISM from the ra- 
dio jet is likely to inflate large cavities (e.g. Mathews & 
Brighenti 2008), which are not seen in NGC 4649 (Humphrey 
et al. 2008). Still, it remains plausible that there could be 
non- negligible cosmic ray pressure in the vicinity of the 
(weak) jets, which extend for <1.5 kpc (Shurkin et al. 
2008). Intriguingly, our current measurement (Fig 5) is con- 
sistent with f n th being significantly weaker outside the cen- 
tral ~2 kpc. In fact, the fit is formally indistinguishable 
if the nonthermal pressure is set to zero at these scales 
(x 2 /dof=19.6/16, as compared to 18.7/13). If cosmic ray 
injection is primarily responsible for the nonthermal sup- 
port in NGC 4649 and similar objects, we would expect to 
see a strong correlation between the radii requiring f n th^ 
and the extent of the jets. This is potentially testable as the 
jet morphology differs widely in nearby galaxies that are 
accessible for this kind of joint X-ray-optical analysis. 

5.2.4 Magnetic Fields 

Observationally, there are known to be > /iG magnetic fields 
in the centres of some early- type galaxies (Vallee 2011), 
which may provide an additional source of nonthermal pres- 
sure. For an early-type galaxy, a plausible mechanism for 
maintaining such a field would be the turbulent dynamo 
effect (e.g. Mathews & Brighenti 1997; Brandenburg & Sub- 
ramanian 2005). Assuming that the entire nonthermal pres- 
sure is provided by this mechanism, and assuming equipar- 
tition between the magnetic and turbulent energy densities 
(which should maintain a disordered field), we computed 
the B-field required in NGC 4649, using P mag = B 2 /8n, 
where P mag is the magnetic pressure and B is the B- 
field in cgs units. The implied field strength profile, shown 
in Fig 6 (left panel), falls steeply with radius (roughly 
B = 24{R/kpc)-°- 8 f j,G) and requires turbulent velocities, 
vturt — 230 km s _1 . The radial B-field dependence is simi- 
lar to what is expected from simple galaxy formation models 
(Mathews & Brighenti 1997; Beck et al. 2012), and compara- 
ble fields have been inferred at small scales in a few galaxies 
with embedded radio jets (Vallee 2011), albeit such strong 
fields are not ubiquitous. 

In the right hand panel of Fig 6, we show how 
the implied B-field varies with the electron density. The 
data were well fitted (x 2 /dof= 1.1/4) with a model B = 
B (n e /0.1cm~ 3 ) 7 , where B = 39 ± 6^G and 7 = 0.59 ± 
0.09. We note that this is similar to the predictions of 
magnetic field generation in a (nonrotating) cooling flow 
predicted by Mathews & Brighenti (1997), albeit the field 
strength is much higher (Fig 6). Faraday rotation measure 
constraints in some galaxy clusters imply 7 ranging from 
~0. 5-1.0 (e.g. Dolag et al. 2001; Guidetti et al. 2008), simi- 
lar to our inferred value for NGC 4649. 

5.3 Measuring gas motions with Astro-H 

If turbulence or bulk gas motions are primarily responsi- 
ble for the nonthermal support in NGC 4649, we should see 
associated Doppler broadening or centroid shifts of the X- 
ray emission lines. Given the spatially extended nature of 



T 




Figure 7. Portions of two simulated Astro-H SXS spectra, show- 
ing a Si XIII line. The spectra were accumulated from the entire 
SXS field of view in two different pointings, which were centred at 
different positions along the major axis of NGC 4649, symmetri- 
cally about the core of the galaxy. The simulations assumed that 
nonthermal support was entirely due to gas rotation; the Doppler 
shifting of the X-ray lines produces a clear shift in the line cen- 
troid between the two spectra (shown as crosses and diamonds). 
For illustration purposes, the simulated spectra shown have 1 Ms 
exposures, as opposed to the 200 ks exposures discussed in the 
text. 

the hot gas, no currently in-orbit instrument would allow 
these modest shifts to be measured definitively. However, 
the nondispersive, high resolution spectroscopy enabled by 
microcalorimeters scheduled to fly on future missions is ide- 
ally suited to this task. 

Of particular interest is the Soft X-ray Spectrometer 
(SXS) aboard the Astro-H observatory (Takahashi et al. 
2010). To explore whether it will be possible to measure 
the predicted gas motions in NGC 4649 with this instru- 
ment, we simulated two deep (~200 ks) exposures, each 
centred along the major axis at a distance ~1.5' from the 
galaxy core. Spectra were extracted from the entire 3x3' 
field of view in each case. If the gas is rotating, one spec- 
trum should be redshifted, and one blueshifted. To simu- 
late the spectra, we modified the Monte Carlo approach 
outlined in § 3.3. At a series of grid positions (R,6>,z), we 
generated photons, assuming a thermally broadened APEC 
plasma with a given density, temperature, abundance, line- 
of-sight velocity and (turbulent) Doppler line broadening. 
These parameters were derived from the best-fitting mod- 
els, assuming spherical symmetry. Assuming isotropy, tur- 
bulence should produce Gaussian broadening of the line, 
/ oc exp (— (E — Eo) 2 /AE 2 ^ , where Eo is the line centre, 

and AE/E = y / 2/3v tU rb/c = b/c, where b is the Doppler 
b-parameter. 

To incorporate the sensitivity of the SXS, we 
used the standard "baseline" on-axis Astro-H re- 
sponses, ah_sxs_7ev_basef ilt_20090216 . rmf and 
sxt-s_100208_ts02um_of _intallpxl . arf 5 . Photons were 
then projected onto the sky, folding in the nominal on-axis 

5 http:/ /astro-h. isas.jaxa.jp/researchers/sim/response. html 
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point spread function (the ' 'ah_sxt_psfmodel_20090217' ' 
model) 6 . Additional background photons were generated, 
corresponding to emission from low mass X-ray binaries (a 
7.3 keV bremsstrahlung model distributed like the stellar 
light and normalized according to the measurement of 
Humphrey & Buote 2008), the cosmic X-ray background (a 
powerlaw with r=1.41 and normalization given by De Luca 
& Molendi 2004), the Galactic foreground (a 0.07 keV and 
a 0.20 keV APEC plasma component, e.g. Humphrey et al. 
2011) and emission from the Virgo cluster ICM (a 2.5 keV 
APEC component; Humphrey et al. 2008). For simplicity, 
we ignored the (featureless) instrumental background which 
is at least an order of magnitude fainter than the source 
emission below ~3 keV. 

Since there was a range of gas temperatures pro- 
jected into the field of view, we modelled the background- 
subtracted spectrum using two APEC models with (the 
same) Doppler broadening and tied abundances and red- 
shifts, plus a bremsstrahlung model to account for LMXBs. 
We found this model adequately fitted the spectra. If gas ro- 
tation supplies the nonthermal support, we would expect a 
peak ~500km s _1 velocity gradient between the two fields, 
which corresponds to ~4 eV at 2.2 keV. In practice, the 
large apertures used caused the line profiles to be smeared 
out. Nevertheless, as shown in Fig 7, the line centroid shift 
should be detectable; for 200 ks exposures, we expect to 
detect the gas rotation at ~4-<r. Conversely, if turbulence 
dominates the gas dynamics, we would expect ~200km s _1 
broadening of the lines, and little velocity gradient. Simu- 
lating corresponding spectra, we found the line broadening 
was required; we constrained the line of sight rms velocity 
to 180 ± 40km s _1 . Therefore, if gas motions are responsible 
for the nonthermal support in NGC4&49, we expect to be 
able to measure them with Astro-H. 
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APPENDIX A: COMPARISON WITH X-RAY 
ANALYSIS OF DAS ET AL. 

Das et al. (2010) introduced a minimally parametric method 
for recovering the mass profiles of systems in which the gas 
is approximately hydrostatic, by inverting the temperature 
and density profile data. Their approach is a variant of the 
"smoothed inversion" technique (for a critical review of mass 
modelling techniques, see Buote & Humphrey 2012a), in 
which, rather than fitting smooth models to the density and 
temperature profile, they instead apply a smoothing prior, 
controlled by the parameter A, which is calibrated against 
realistic data. They assumed that A is the same (and fixed) 
for all observables (temperature, density, rotation velocity) 
and all systems. The advantage of this method is that it 

6 http: / /astro-h. isas. jaxa. j p /researchers / sim / GSFC_mirror . html 



makes few a priori assumptions about the shape of the mass 
or temperature profiles. The disadvantage is that it assumes 
that all of the profiles are smooth in a particular manner, 
which is the same for all profiles 7 . 

In Fig 1, we compare the v c profiles obtained by our 
method and by Das et al. (2010). In general, the agreement 
is good over the range ~3-9 kpc, but the mass found by 
Das et al. (2010) is higher outside of this range. In Fig Al, 
we compare the temperature and density data which were 
fitted in our analysis, and those used by Das et al.. It is 
immediately clear that there is good agreement between all 
three data-sets between ~3-9 kpc, where the v c profiles also 
agree, but discrepancies arise outside this range. We con- 
clude that the differences between the mass inferred from 
the hydrostatic analysis of Das et al. (2010) and Humphrey 
et al. (2008) are most probably a consequence of the dif- 
ferent temperature and density profiles, rather than intrinsic 
differences in the mass fitting techniques. 

At the smallest scales, there is a large difference be- 
tween the XMM and Chandra data-points, which proba- 
bly reflects the larger PSF of XMM (90% encircled energy 
radius~30"=2.2 kpc) compromizing the inner data by spec- 
tral mixing between adjacent annuli. This offset may explain 
the formally poor reduced \ 2 (—6) obtained by Das et al. in 
their fits. We note that such large systematic errors will dis- 
tort the x 2 topology, even if the best-fitting model is correct, 
compromizing the error-bar calculation. 

The remaining differences with our Chandra profiles (in 
particular the lower temperature we found at large scales) 
probably reflect different analysis choices made by Churazov 
et al. (2010a), who derived the profiles used by Das et al.. 
Specifically, in their analysis the ISM abundance was fixed 
to a constant value (ZF e =0.5), and Solar abundance ratios 
were used for all species, which is not formally correct for 
NGC 4649 (Humphrey & Buote 2006; Humphrey et al. 2006; 
B09). Furthermore, they ignored the contribution of unre- 
solved LMXBs (which can affect the inferred density) and 
projected emission from the Virgo ICM (which can affect 
the temperature) in their fits, while their correction for pro- 
jected emission from gas beyond the outermost shells used 
in the deprojection did not account for the truncation of 
the halo by the Virgo ICM. In our modelling, we explicitly 
accounted for these effects (Humphrey et al. 2008). Never- 
theless, even with the different analysis choices, we see that 
the recovered mass profile is not dramatically affected over 
most of the interesting radial range (Fig 1). 
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